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Since its discovery in 1995, the top quark properties were extensively studied at the Fermilab 
Tevatron. In the last 15 years, many observables concerning top-quark physics were measured 
with remarkable accuracy. Among others, the tt total cross section. Off, was measured with an 
accuracy of AOtt/Oft ~ 9%. Data coming from the LHC are expected to improve significantly the 
measurement of several observables related to the top quark. Already within a couple of years 
of data taking in the low-luminosity low-energy phase ~ lOOpb^^/year at 7TeV of center of 
mass energy), tens of thousands tt events before selection will be available. Consequently, already 
in this first phase, the accuracy of the cross section measurement is supposed to match the one 
reached to date at the Tevatron. In the high-luminosity phase ( ^ ~ lOOfb^^/year at 14TeV of 
center of mass energy) it will be possible to reach an accuracy AGn/On ~ 5% |[T]]. This accuracy 
in the experimental measurements motivates theorists to refine the existing predictions, both for 
the total top-quark pair production cross section and for the related differential distributions. In 
the following, we will briefly outline the research program aiming to the full calculation of the 
next-to-next-to-leading order (NNLO) corrections to top-pair production cross section. 

The full NLO QCD corrections to the total cross section were calculated in |§] in the case of 
"stable" on-shell top-quarks. In ^ several differential distributions were calculated at the same 
accuracy level. In [Q] the NLO corrections to the top-pair production were evaluated by taking into 
account the top-quark decay in narrow-width approximation. The resummation of soft-gluon en- 
hanced terms near the production threshold is implemented at the leading [||], next-to-leading [^], 
and next-to-next-to-leading logarithmic accuracy. Approximate NNLO formulas for the total 
cross section were recently obtained by several groups [^. The calculation of the corrections to the 
partonic cross section beyond leading order can be split in the calculation of the real corrections, in 
which the final state includes extra partons in addition to the top-quark pair, and in the calculation 
of virtual (loop) corrections to the partonic processes already present at the tree level. A complete 
calculation of the NNLO corrections requires the knowledge of the two-loop matrix elements for 
the processes qq — )• tt (quark-annihilation channel) and gg tt (gluon-fusion channel), as well as 
the 2 — >■ 3 matrix elements at the one-loop level and the 2 — 4 matrix elements at the tree-level |^]. 
Moreover, in order to be able to deal with IR singularities, a NNLO subtraction method has to be 
implemented [llO|]. 

From the technical point of view, the calculation of the two-loop (virtual) corrections is par- 
ticularly challenging. The squared matrix element for the processes q{pi ) + q{p2) — ^ l{P3) + l{P4) 
and g{pi ) +g{p2) — > t{P3) + t{P4)^ summed over spin and color, can be expanded in powers of the 
strong coupling constant as as follows: 



'{s,t,m,e) = 1671^ a| 



(1) 



where = for / = 1,2 and pj = —mj for / = 3,4. The Mandelstam variables are defined in 
the usual way: s = — {pi + P2)^, t = —{p\ — p^)'^, u = —{p\ — p^xf . Conservation of momentum 
implies that s + t + u = 2mj. 

The tree-level term £/o in the r. h. s. of Eq. ([T]) is well known in both production channels. The 
I0'{as) term arises from the interference of one-loop diagrams with the tree-level amplitude. 
The ff{(xj) term consists of two parts: the interference of two-loop diagrams with the Bom 
amplitude and the interference of one-loop diagrams among themselves, £^2 = i*^'^^""* 
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^(1x1) 



.(2x0) . . , 

2 , origmat- 



The latter term, , was calculated for both channels in [|lip. The first term, 

ing from the two-loop diagrams, can be decomposed according to the color and flavor structures as 
follows: 



^^(2x0)?^-^ 



^i^m^B'i^^ ^ +N,NiDf+ ^Ef+NMDf+ ^Ef+Nfpf 



+NtN,Fl^+NlFf 



, (2) 



where Nh and A'^^ are the number of heavy-quark flavors (in our case, only the top quark is consid- 
ered heavy) and light-quark flavors, respectively. The coefficients A,B,.. .,///, in both channels are 
functions of s, t, and m,, as well as of the dimensional regulator e. These quantities were calcu- 
lated in 1 12] in the approximation s, |m| ^ mj. However, the results in [ p^ ] are not sufficient to 
obtain accurate NNLO predictions, since a large fraction of the events are characterised by values 
of the partonic center of mass energy which do not satisfy the ultra-relativistic limit. The com- 
plete top mass dependence of the color coefficients A, B, ... is required. A numerical calculation of 
, exact in s, t, and mt, was presented in [|l^]. Analytic expressions for all of the IR poles 



(2x0) qq 



m 



^(2x0)?^ 



and 



^(2x0)gg 



are also available [ p3| ] ; they were calculated by employing the expres- 
sion of the IR divergences in a generic two-loop QCD amplitude with both massive and massless 
particles, derived in [[l4|]. 

The analytic calculation of the coefficients A^^, Df ,...,F^'^ which appear in the quark anni- 
hilation channel, as well as of the coefficient A^^ for the gluon fusion channel, was carried out in 



[16]. It must be observed that the coefficients with a I or h subscript receive contributions only 
from diagrams involving a closed light or heavy fermion loop, respectively. The leading color co- 
efficients in the two production channels, A^^ and A^^, involve planai" diagrams only. The results 
reported in [16] were obtained by employing the Laporta Algorithm [p^, implemented in the C++ 
code REDUZE 2 [jl8|, |l^], for the reduction to the Master Integrals ^ [T^. Subsequently, the 



master integrals were evaluated by means of the differential equation method []21p. The analytic 
expression of the master integrals can be written in terms of one- and two-dimensional harmonic 



polylogarithms [|22|]. The analytic results were evaluated numerically by employing codes which 
make use of a GiNaC package for the evaluation of generalized polylogarithms [E3]. In [16] the 



analytic results were also expanded in the s » mj limit, in order to reproduce the results already 
obtained in [ 12]. Starting from the exact results of [p^], it was also possible to obtain new analytic 
formulas which are valid in the production threshold limit s — )■ 4mj. 

A complete numerical result for the two-loop corrections in the gluon fusion channel is still 
missing, and to date only the coefficient A^^ was evaluated. Among the remaining color coefficients 
in Eq. (^, Ef^, Ff^, Gf, Hf^, Hff , Hfjf, if^, if^^ , /f^f can be calculated using the same technique 
already employed in [p^; their evaluation is in progress [Q]. The remaining color coefficients both 
in the quark-antiquark channel {B^^ and C^^) and in the gluon fusion channel (S^^, C^^, D^^, Ef^^ , 
F^^, and Gf^) contain either crossed box topologies, or complicated massive sub-topologies. In the 
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first case, the calculation of the color coefficients in a closed functional form using the differential 
equation method is very difficult, because of the large number of master integrals that occurs for 
some specific topologies. In the second case, the problems arise from the fact that it is known that 
already the sunrise diagram with three virtual propagators of mass m and an external momentum p 
such that 7^ can be expressed analytically only in terms of elliptic integrals [^]. This three- 
denominator diagram appears as a sub-topology in the color coefficients with subscript h. A viable 
solution for both issues is the semi-numerical approach adopted, for instance, for the two-loop 
equal-mass sunrise diagram [26|. 

To conclude, the calculation of the two-loop corrections to the top-quark pair production is an 
essential step needed for the full evaluation of the NNLO corrections to the top-quark production 
cross section and differential distributions. In the last few years several results were obtained for 
the quark annihilation channel. The calculation of the two-loop corrections in the gluon fusion 
channel is technically more complicated because of the larger number of diagrams involved and 
because the functional basis of the harmonic polylogarithms is known to be insufficient to obtain 
a full analytic result. However, the calculation of a large number of the color coefficients can be 
carried out with available methods and is in progress. The calculation of the color coefficients 
for which standard analytic techniques are insufficient can in principle be carried out with semi- 
analytic methods already applied to other related problems. 
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